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ABSTRACT 

We use the (^-function regularization and an integral representation of the 
complex power of a pseudo differential operator, to give an unambiguous defi- 
nition of the determinant of the Dirac operator. We bring this definition to a 
workable form by making use of an asymmetric Wigner representation. The ex- 
pression so obtained is amenable to several treatments of which we consider in 
detail two, the inverse mass expansion and the gradient expansion, with concrete 
examples. We obtain explicit closed expressions for the corresponding Seeley- 
DeWitt coefficients to all orders. The determinant is shown to be vector gauge 
invariant and to posses the correct axial and scale anomalies. The main virtue of 
our approach is that it is conceptually simple and systematic and can be extended 
naturally to more general problems (bosonic operators, gravitational fields, etc). 
In particular, it avoids defining the real and imaginary parts of the effective action 
separately. In addition, it does not reduce the problem to a bosonic one to apply 
heat kernel nor performs further analytical rotations of the fields to make the 
Dirac operator Hermitian. We illustrate the flexibility of the method by studying 
some interesting cases. 
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1. Introduction 

The success of the Lagrangian formulation of relativistic field theories is due to the 
easy implementation of symmetries and in particular of Lorentz invariance. The same is 
true for the associated covariant quantization procedure, the Feynman functional integral 
approach. In a typical setting, as will be ours in this work, one considers Dirac fermions in 
D space-time dimensions in the presence of arbitrary non gravitational bosonic background 
fields coupled to Lorentz and to internal symmetry indices of the fermions. Generally 
speaking, the background fields can be treated as external. This is no restriction since 
they can be quantized introducing the corresponding functional integration over them. 
The effective action of the system is then obtained by integrating out the fermion fields. 
Formally the Grassmann integral gives the determinant of the Dirac operator, DetD, and 
the effective action is just its logarithm. As a consequence this determinant plays an 
important role in the functional integral formulation [1]. 

In the relativistic case the naive determinant is ultraviolet divergent and one has to 
define a renormalized determinant by introducing counterterms or some other equivalent 
technique. In this process some classical symmetries can be lost. Of course, all this parallels 
the diagrammatic approach where the determinant is represented by the one fermion loop 
graphs. 

There has been a number of ways to address the problem of defining a finite de- 
terminant or equivalently a finite effective action [2]. Most of the methods developed in 
the literature try to reduce the problem to a bosonic one. This is because second order 
differential operators are algebraically simpler and better studied in this context. In addi- 
tion, the important ultraviolet problem can be treated within an inverse mass expansion 
(i.e. a simultaneous weak and smooth field expansion) with the heat kernel technique, 
and this applies to second order definite positive Hermitian differential operators. The 
reduction is usually achieved by considering or DD^ and D^D [3-8]. In the first case 
a further analytical rotation of the fields (besides the Euclidean rotation) is assumed to 
make D antihermitian. This requires to extend the internal symmetry group too and it 
is known that the rotation back of concrete subgroups can be ambiguous [2]. Further- 
more, although looks like a bosonic theory, some of the efficient methods developed 
to go beyond the inverse mass expansion [9], may not be straightforwardly applied to 
fermions. This is because such methods assume the unrestricted validity of the formal 
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relation Det{AB) = Det (^)Det (S), which does not hold under regularization. In the 
bosonic case the use of this formal relation just redefines the counterterms, but for fermions 
it makes the determination of the chiral anomaly ambiguous. Other approaches use for- 
mal relations to define independently the real and imaginary parts of the effective action, 
using D and D''' and are rather ad hoc [8,10]. Similar approaches used to define an ad 
hoc Jacobian of the functional measure under chiral transformations are known to lead to 
results inconsistent with the Bardeen anomaly [11]. 

A diff'erent approach is that of Leutwyler [12] and Ball [2]. It uses a formal definition 
of the variation of the effective action in Euclidean space which is regularized in a chirally 
invariant way using D^D in a proper time representation. For the real part of the action 
this defines a true variation. The regularized variation of the imaginary part satisfies the 
integrability condition only after adding suitable polynomial counterterms which introduce 
the chiral anomaly. This approach is both mathematically impeachable and also computa- 
tionally convenient in heat kernel-like expansions. However it is rather sophisticated and 
relies heavily on the previous knowledge of the subtleties involved in Det D obtained after 
years of deeply original insights [13-20]. For instance the operator D^^, which in principle 
is unrelated to the problem and formally cancels in the definition, is needed in the con- 
struction. It is not clear to us how this approach can be properly extended to more general 
theories, such as non local theories, string theories, and so on. 

We think that it can be of interest to pursue an alternative statement of the problem 
where the definition of the determinant is given at the very beginning and then every 
other quantity can be defined and calculated unambiguously without the need of new 
prescriptions. A good example of the latter is provided by the functional Jacobian under 
variations of D. Although such an object is not needed in this kind of approach, it is 
perfectly well defined and can be computed if desired [21]. The only remaining freedom is 
of course the addition of counterterms, polynomials of degree D in the number of derivatives 
plus external fields, which allows to reproduce any other renormalization prescription or 
enforce particular symmetries. A suitable definition, given long ago [22], is the C- function 
regularization which has the advantage of preserving automatically a large class of the 
classical symmetries, namely those which are implemented by similarity transformations 
of D. This includes relativistic invariance, vector gauge invariance and so on. Other 
virtue is that it is well grounded mathematically [23]. Using a typical Cauchy integral 
representation, the ^-function can be related to the resolvent of D which is a more tractable 
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object and is suitable to perform systematic expansions. Finally, a Wigner representation 
[24,25] allows us to treat properly the ultraviolet divergences appearing in the diagonal 
elements of the resolvent. In this way a well defined workable form for the effective action 
and other quantities is obtained. Although by no means is it intended in this work to 
achieve strict mathematical rigor, the troublesome ultraviolet sector is treated carefully, 
and only in the infrared we proceed formally by assuming that D is well behaved in this 
sector, that is, effective boundaryless boundary conditions and no zero modes. 

Because the subject has been extensively studied in the past, we do not intend to 
present truly new theorems, rather our emphasis is on introducing a conceptually simple 
scheme to make the subject more easily graspable with more systematic, and sometimes 
simpler, proofs of known results and also with an eye put on the generalization to worse 
known systems where our intuition is less developed. 

The paper is organized as follows. In Section 2, we present our conventions and set 
up the basic formalism for the ^-function regularized determinant of the Dirac operator, 
as well as the corresponding consistent currents. Actually, we find that the effective ac- 
tion can be reconstructed from the current without loss of information in this particular 
regularization. Section 3 deals with one special definition of the Wigner transformation, 
particularly adapted to the Dirac operator, and rather convenient from a computational 
point of view. In Section 4 we consider an inverse mass expansion for the Dirac opera- 
tor and among other things, we obtain a direct, i.e. non recursive, determination of the 
corresponding Seeley-DeWitt coefficients. These can be profitably used to write down in- 
verse mass expansions both for the effective action and the effective currents. A further 
application of the previous results can be found in Section 5, when computing chiral and 
scale anomalies within the ^-function regularization. Moreover, we establish the general 
form of the counterterms needed to bring the chiral anomaly to its minimal (Bardeen) 
form. To do so the Wigner transformation method turns out to be very useful. In Section 
6, some aspects of the so-called gradient or derivative expansion are studied within the 
present formalism, i.e. Wigner transformation technique. We consider for instance a direct 
calculation of the gauged Wess-Zumino-Witten term from the Dirac determinant. We do 
so without using better known methods, like e.g. trial and error, differential geometry or 
integration of the chiral anomaly. Along similar lines, we study fermionic currents and 
their relation to effective actions in two less dimensions, again with the help of the Wigner 
transformation method. Finally, parity anomalies in odd dimensions are revisited in Sec- 
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tion 7, where the gauge invariant Chern-Simons action is obtained. In Appendix A, we 
estabhsh a exphcit closed formula for the Seeley-DeWitt coeflBcients of the Dirac operator 
to all even orders. In Appendix B, the connection between the well-known Heat Kernel 
expansion for the squared Dirac operator, D^, and our inverse mass expansion for the 
Dirac operator, D, is developed in the even-dimensional case. Appendices C and D collect 
explicit formulas for the chiral anomaly and the counterterms for the most general Dirac 
operator in four dimensions. 

2. (^-function regular izat ion 

Let D be the Dirac operator in D dimensional Euclidean space for a fermion in the 
presence of arbitrary (non gravitational) external fields 

D = z^+Y(a;) (2.1) 

Here Y(a;) is a matrix in the internal degrees of freedom, i.e. spinor and flavor, but does 
not contain derivative operators. Our conventions regarding gamma matrices are as follows 

ll = -7m , {7m> li^} = -2<^M^ ^2 2) 

75 = -*^^^7o7i • • • iD-i , eoi...D-i = +1 

Whenever needed we will assume the standard hermiticity for the external fields, that 
is, such that if they transform covariantly under the Wick rotation, 70D is Hermitian in 
Minkowski space. This implies that the Euclidean effective action is real in the pseudopar- 
ity even sector (containing no Levi-Civita pseudotensor) and imaginary in the pseudopar- 
ity odd one (containing a Levi-Civita pseudotensor). We will often use the object D itself 
rather than Y{x) because this produces more compact formulas. Another important point 
is that D transforms homogeneously under the classical symmetry transformations 

D ^ D" = l^aDl^j;^ 

i;^{x) = Qii;{x') (2.3) 
^"(x) = V;(xOl^2"' 

whereas Y, in general, transforms inhomogeneously. This transformation corresponds to 
a classical symmetry if it maintains the structure (2.1), i.e. = i ^ +Y^(a;), but 

otherwise fii, and Q.2 can depend on x and contain derivative operators. In this paper 
we will consider explicitly two classical symmetries, namely chiral gauge rotations and 
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scale transformations. Notice that the form of the Dirac operator does not include general 
coordinate transformations as classical symmetries. This would require to extend the 
Dirac operator by properly coupling gravitational fields, and subsequent generalization of 
our computational procedure to curved space-time. Such a study will not be undertaken 
here and is left for future research. 

For definiteness we can think of D as admitting a complete set of left and right 
eigenvectors 

Dcpnix) = Xn(l>n{x) , Cn{x)'D = XnU{x) (2.4) 

which can be normalized (in a box) so that (Cnl^fc) = ^nk- More generally D can have 
a general Jordan form, i.e. completeness of the eigenvectors will not be required. The 
Euclidean partition function is a functional of the external fields Y{x) given by 
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Vi(;Vi(;ex-p!^- J d^xij{x)'Dij{x)^ (2.5) 



and formally Z = Det D. Hence the Euclidean effective action is formally 

= - log Z = -TV log D = - J] log(A„) (2.6) 

n 

The sum in (2.6) is ultraviolet divergent and must be regularized. To do so we shall 
adopt the ^-function regularization prescription [1,22]. We shall assume that D has no 
eigenvalues in some neighborhood of zero, otherwise some infrared regularization would be 
needed too. If there is only a finite number of zero modes, isolated from the rest of the 
spectrum, one can define a restricted determinant excluding the zero modes. 

Let us consider the pseudo differential operator (D///)*, where is a scale introduced 
for dimensional bookkeeping and the complex number s is the regulator. If D admits a 
complete set of eigenvectors, this operator is characterized by its eigenvectors and eigen- 
values {(f)n, (An/A*)*}- III case a convenient representation is given by [23] 

/Dy _ r dz /zy 1 

where F is a path that starts at infinity, follows a ray of minimal growth (i.e. no eigenvalue 
of D lies on it), encircles the origin clockwise and goes back to infinity along the ray. No 
eigenvalues are encircled by F. If the set of eigenvalues is bounded, F can be deformed to 
an anticlockwise oriented closed path containing the eigenvalues of D but excluding the 
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origin. This construction assumes that no eigenvalue hes in a neighborhood of the origin 
and hence that D is non singular. This representation is meaningful for Re s < 0, and can 
be analytically continued to the complex plane of s. With this prescription, (a;|D*|?/) is an 
entire function of s for x ^ y and has simple poles at s = —1, —2, . . . , —D, foTX = y [23]. 
(Here \x) is an eigenvector of the multiplicative operator x^ in D-dimensional Euclidean 
space. In what follows we will not use a different notation for a function of and 
the corresponding multiplicative operator, and similarly we will identify the momentum 
operator and id^). In particular D'^ is analytic at s = and this allows for the definition 
of a regularized logarithm of D 

log(H) d.Dy 

\ /I / as \ ji / s=o 

and of the regularized determinant of D [1,22] 

W(D) = -logDet (^) = -Tr log(^) = -AJ d°.tr(x|(^) (2.9) 

The determinant so defined is completely finite but depends on the arbitrary scale /i 
introduced by the regularization. 

Using previous formulas it is easy to prove the following identity satisfied by the action 
in this regularization 

d ris + 1) f dz d 



ds r(s + n) Jy 27rz oz 

In fact in the r.h.s. we can use any other regularization because all of them differ by a 
polynomial in z and its contribution cancels in the integral. On the other hand using a 
value n > D, d^WCD — z) is ultraviolet finite and independent of the regularization and 
can be used to reconstruct the action. 

It is also of interest to obtain expressions for the variation of the action under a 
generic infinitesimal transformation of the fields in the Dirac operator. Let 5xD = X be 
such a variation of D, where the only restriction is that X is a multiplicative operator, 
that is, without derivatives. Hence it corresponds to an infinitesimal classical symmetry 
transformation. Correspondingly 

<5xD = X, 5xW^= / d^xtr(XJ), (2-11) 

We will refer to J as the current associated to W, although usually this name is reserved 
for the variation under gauge fields and otherwise J is called a density. The current so 
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defined satisfies consistency conditions [17] which reflect that it is a true variation and 
hence it is known as the consistent current. Other definitions of the current are sometimes 
more convenient, in particular the so caUed (chiral) covariant current [26], which wiU be 
considered in Section 5. They differ by a polynomial in the fields and derivatives. 

An explicit expression for the current can be obtained by using the technique intro- 
duced in [21,27] 



d 

as jj^ 



S-1 



(2.12) 

s=0 



s=o ds 

We have used the cyclicity of the trace in presence of the regulator s. Because X is local, 
we have 

d D*"^ 

J(x; D) = - — s{x\ \x) (2.13) 

The current can be used to recover the action by applying Eq. (2.10) for n = 1, 

W{-D) = -j d^'xtv^ j^^yj{x;n - z)\s=o (2.14) 

Again, we can use any other regularization for the current in the r.h.s., moreover there is 
no need to use the consistent current, it can be for instance the covariant current, and the 
z integral will take care of introducing the proper chiral anomaly into the action. 

The ^-function regularized action enjoys all the classical symmetries which are also 
symmetries of the quantum Dirac equation i.e. the eigenvalue equation = Xn(pn- 

This includes in particular vector gauge invariance, D — > Q{x)JDQ~^{x), with Q{x) a 
matrix valued function acting on flavor space but not in Dirac space. On the other hand 
scale and axial transformations are symmetries only of the classical equation = 0, 
and the corresponding currents are anomalous. Note that the classical equation does not 
imply a zero mode in the quantum equation due to the different boundary conditions 
of both equations. Actually the solutions of the classical and the quantum equations 
are normalizable in a spatial box and space-time box respectively. The scale and axial 
anomalies will be considered in Section 5. 

3. Wigner transformation 

Seeley's representation (2.7) requires to invert the operator T) — z. This can be con- 
veniently accomplished by means of an asymmetric version of the Wigner representation 
[24,25]. For any operator A, let 

A{x,p) = J d''ye^y^x\A\x-y)=^-^^ (3.1) 



be its (asymmetric) Wigner representation. \p) is the momentum eigenstate with {x\p) = 
g-iap When ^ is a pseudo differential operator [28], the quantity so defined is closely 
related to the so called symbol of A. For a pseudo differential operator A of complex order 
m, A{x,p) is of the form '^^Qam-j{x,p) where ak{x,p) are homogeneous functions of 
degree k in p. In particular D* is of order s [23] . 
From this definition 

-^^A{x,p), TrA = J -^^-^trA{x,p), (3.2) 

where tr acts on internal and Dirac spinor degrees of freedom only, and the product of two 
operators satisfy the following formula 

{AB){x,p) = exp(za^ • d^)A{x,p)B{x,p) , (3.3) 

where dp acts only on the p-dependence in A{x,p) and on the x-dependence in B{x,p). 
Let the propagator or resolvent of D, be 

G{z) = (D - z)-' (3.4) 

and G{x,p]z) its Wigner representation. A convenient expression for G{x,p;z) can be 
obtained using the trick of ref. [9]. Recalling Xi^\x) = x^lx) 

(x\G(z)\p) = (a;|e-*'=*e+^'=*(D - z)-'^\p) 

(3.5) 

= e-*'=^(a;|(^ +D - z)~^ e+'''^ \p) = {x\k){x\{^ +D - z)-^\p - k) 

where we have made use of e^'^^^ id ^e''^^^ = id^ + k^. Choosing k = p we obtain the 
following compact expression for the Wigner representation of the propagator 

G{x,p;z) = {x\{^+n- z)-^\0) (3.6) 

where |0) is the state of zero momentum, (a:|0) = 1. In practice this implies that 
derivates every x dependence at is right, until it annihilates |0). At this point we have 
separated from the standard approach, which uses the product formula (3.3) to set up 
a recurrence relation to compute the symbol of the resolvent [21]. Our method is more 
efficient for it computes directly, that is non recursively, each of the terms. 

Several expansions can be devised to compute G{x,p;z). Two of them will be con- 
sidered in the next sections. Let us point out that the definition given for A{x,p) is not 



9 



gauge covariant because \p) is not. It would be very interesting to have a gauge invariant 
version of the Wigner representation. In fact the problem of finding non local covariant 
expansions, as required for instance in massless theories, is still open [2, p. 128]. In what 
follows we will consider only local objects of the form (a;|/(D)|a;) as given by the formula 

We will assume that the function f{z) is sufficiently convergent at infinity or else that it 
can be obtained as a suitable analytical extrapolation from a parametric family f{z, s) in 
the variable s. In either case the integration over z should be performed in the first place, 
to yield the Wigner representation of the operator /(D). Afterwards, the p integration 
is carried out, corresponding to take the diagonal matrix elements of (y|/(D)|a;), hence 
restoring gauge covariance. 

Because the subject of this paper has been considered extensively in the past [2] , and 
to some extent it has been taken up by mathematicians, we must make some comment on 
the validity of equations such as (3.2) from a more rigorous point of view. This equation 
is based on an assumption, namely that there exists an operator £^ such that it satis- 
fies Heisenberg commutation relations [id^^Xy] = i5^y. Of course this is not the case for 
manifolds without a global chart and in particular for the nice compact manifolds without 
boundary, usually considered in the literature. For instance if the previous formalism is 
applied to compute Tr (D^) with D = id^ + m on the interval [0, T] compactified to a 
circle, we will find that the sum on the eigenvalues of the exact formula has been replaced 
by an integral on a continuum momentum label, thus introducing an approximation which 
becomes exact as T ^ oo. On the other hand there is no problem for if the compactifi- 
cation comes from the nature of the external fields in the infinity. For instance an equation 
similar to (3.2) for the one dimensional second order differential operator A = —d^ + 
would still be exact, including detailed information on the discrete spectrum of A. How- 
ever such a formula would not be very useful unless one makes some expansion, typically 
a gradient expansion (to be considered in Section 6), and this kind of expansion substi- 
tutes the discrete spectrum by a smoothed continuous density of states. This is a common 
feature of any asymptotic expansion such as heat kernel [29] or the Wigner-Kirkwood [30] 
expansions. 

4. Inverse mass expansion for Dirac operators 

4.1 Seeley-DeWitt coefRcients for Dirac operators 
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In this section we will consider a Dirac operator D = Dq + m, with m a mass term 
and Do of the general form (2.1). This will allow us to obtain an expansion for the effective 
action in powers of Dq or equivalently in inverse powers of m. This expansion is closely 
related in spirit to the Seeley-DeWitt expansion for the heat kernel of positive definite 
second order operators. Both are asymptotic expansions in the number of external fields 
and their derivatives. We will consider the quantity {xKY)/ n)^\x) from which the effective 
action is easily obtained, Eq. (2.9). From Eq. (3.6) it is straightforward to derive the 
following expansion 

where recall that 9^ inside Dq annihilates |0). 

The series will generally be asymptotic, giving only the analytical part of {x\{]/) +m — 
z + ADo)~"'^|0) as a function of A about A = 0. This expansion is thus relevant for large 
Pn, z and m or equivalently for sufficiently weak and smooth external fields. 

Inserting the series for G{x,p; z) in Eq. (3.7) for f{z) = {z/ jiY^ the following expan- 
sion will be found 



(■r-j\ s °° 1 s+D—N 

-) I-) = E (4,)i./2" Ms){x\0^iT^om (4.2) 



where D is the space-time dimension, cjv(s) is a numerical function of s, N and D, and 
On (Do) cire polynomials of degree N in Dq which can be written in a D independent form 
and are matrices in internal and Dirac spinor space. The precise factorization between 
cn{s) and On is, to some extent, a matter of convenience. 

To obtain explicit expressions for On it is best to work out the numerator in Eq. (4.1) 
using the object = ^{7^, Dq}, and the property [p^A^,|^] = 0. Afterwards, an angular 
average over is done so that the following integral applies: 

/i(fc,a,6,A^,m) = J j^^yp^\z-mr{p^ + {z-mf)-' (4.3) 

where m > 0, F goes along the real negative axis, -D, /c, a, h gZ, L>, a > 0, -D + 2/c > and 
Re{s) < 2h - D - 2k - a - I. Exphcitly 

_ {-l)^ F(iD + A-)F(6- iD-A0F(.s + l)/;r^+i+^+2'''^"-2^ 
^1 " (47r)^/2 V{\D)T{h)T{2h - D - 2k - a)T{s + 2 + D + 2k + a - 2b) ^ ' ^ 
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If m is allowed to be negative an additional factor e(m)^ has to be included, e{x) being 
the sign of x. This circumstance will be relevant in the odd dimensional case (see Section 
7). 

In this way one obtains for cn{s) 



Cn{S) — (^N-D: 



T{s + 1- N + D) 
C -( 1)^ r(iV-ll^) AT-l^ + ll 



(4.5) 



where \_x\ stands for largest integer not exceeding x. This explicit expression for (^k holds 
for N > D and an analytical continuation in D is understood for N < D. Note that 
depends also on D (mod 2). These quantities can be obtained more conveniently from the 
recurrence 

Ck-1 = ^kCk , 

Co = l (evenD), Co = -r(^) (odd L>), (4.6) 
eK = 2 (evenD + K), sk = K {odd D + K) 

In particular (k vanishes for negative K if D is odd. Some useful particular values of the 
coefficients are given in Table 4.1. 

On the other hand, for the operators Oat at lower orders, one obtains 

Oo(D) = 1 

0i(D) = D + A^7^ 

C»2(D) = D2 + A^ 

03(D) = + ^D^A^T^ + ^A^7^D2 + A^DA^ 

+ -(A2a,7,+ A^7^ + A^A^7^) 

e»4(D) = + ^(D^AJ + A^D^A^ + A^D^) 



2 3 
1 

6' 



~l~ « (-^/n-^i^ (A^Aj/) + A^Aj^A^) 



where A^ = |{7^,D}. We want to emphasize that these expressions do not make any 
assumption on D other than Eqs. (2.1) and (2.2). Note also that they have been written 
in a D independent fashion. The expression for the coefficients Ojv(D) for all even orders, 
is given in Appendix A. 
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In the expansion (4.2), what actually appears are the matrix valued functions On{x) = 
{x\On\0), rather than the operators On themselves. Because the regularization is vector 
gauge invariant, we expect that ONix) will be covariant under vector gauge transfor- 
mations of D, that is, under D — > Q{x)'DQ~^ (x) , On{x) will transform as ONix) — > 
il{x)ON{x)fl~^{x), even if |0) itself is not gauge invariant. Let us show this explicitly for 
02- Letting T) = iTj^ +X with = id^ + Vf^{x), the operator O2 can be written as 

(4.8) 

Therefore all the derivatives appear inside commutators and O2 is a purely multiplicative 
operator. In this case taking the matrix element (a:|(92|0) does not break gauge invariance. 
The same thing can be shown for higher orders. A more economical manner of establishing 
the gauge covariance of an expression is the following. Let / be an operator formed 
algebraically out of iV^ and other gauge covariant non derivative operators Xi(a;). Then 
of course / is itself a gauge covariant operator. The matrix valued function f[x) = (a;|/|0) 
will be gauge covariant if / is a multiplicative operator, that is, if all the in / appear as 
covariant derivatives inside commutators. In turn this will be only the case if / is invariant 
under Vf^{x) — > Vf^ix) + a^, with an arbitrary constant c-number, that is if 

5gf = under ^^D ^^X^ = (4.9) 

For the expressions listed above one checks, using 8gA^ = — a^, that in fact 6gON = 0. 
Note that this rule refers only to vector gauge invariance of local objects. For instance 
in the Chern-Simons action j^i J d^x e^;^Q, tr(V^9iyVa — |«V^Vi^Vq,), (to be discussed in 
Section 7), the integrand is only invariant up to a total derivative but the action itself is 
invariant. 

Another remark about gauge invariance is that it follows from formal integration by 
parts over p [9] in Eq. (3.7). In the case of the integral Ii{k,a,b,D,s,m) in Eq. (4.4), 
only the relationship 

/i(/c + l,a-2,6,A8,m) ^ + k 

Ii{k,a,b,D,s,m) b-\D-k-l ^' ' 

is required to form the gauge covariant operators On- Because this ratio is independent 
of s, the same quantities will appear in the expansion of {x\f{D)\x) for other analytic 
functions f{z) whenever the z integrals involved are convergent. This argument is made 
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rigorous in [23] for f{A), A being a zeroth order pseudo differential operator and f{z) 
analytic in a region containing the spectrum of A. 

It is worth noticing that the operators On are related among themselves in a simple 
way. This follows from the observation that D is invariant under 

Sm^o = —Sm, = 5m (4-11) 

Imposing 5m{x\Ti^\x) = in Eq. (4.2), one finds 

5mON = -eN-oON-iSm (4.12) 

where the quantities ck were introduced in Eq. (4.6). Therefore from On one can obtain 
algebraically Om, M < N. It is interesting to notice that the variations Sg (see Eq. 
(4.9)) and also 5^. (see Eq. (4.11)) greatly restrict, but do not completely fix, the actual 
values of the coefficients appearing in the expression for Cjv as given by Eq. (4.7). It turns 
out that the even order coefficients can be deduced entirely as combinatorial factors, the 
odd orders being obtained by taking the variation with respect the mass, as given by Eq. 
(4.12). Further details can be looked up at Appendix A. 

The coefficient cjv(s) has simple poles at s = —1, —2, —3, . . . , N — D for N < D. This 
is in agreement with general theorems for pseudo differential operators [23] . On the other 
hand cn{s) vanishes for s = 0, 1, 2, . . . , AT — D if A?" > D. Thus for non negative integer 
values of s, the inverse mass expansion is exact [23] and is just a polynomial in m. Setting 
m = 0, one finds: 

(xlD'^lx) = j^^CnTin + l)OD+n{x; D), n = 0, 1, 2, . . . (4.13) 

This result is regularization dependent, since a naive evaluation of the matrix element 
would diverge. In fact, the l.h.s. stands for {x\T)'^~^^\x)\s=o through an analytical contin- 
uation in s, and thus it is specific of the ^-function regularization. 
4.2 Effective action 

For the effective action we have, using (2.9) and (4.2), 

W{1D) = f; + /3r^log (^)) m^-^(0^(Do)) 

iv=o ^ V/^// (^4J4) 

= -^^cn+K{0), = -cd+k{0) 
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In the r.h.s. we have used the shorthand notation 

(/(D)) = / d^^tr(x|/(D)|0) = ^^-|^tr(0|/(D)|0) (4.15) 

The symbol ( ) enjoys some of the properties of the trace and in particular the trace cyclic 
property, but only for multiplicative operators, 

{AB) = {BA) if 5gA ^6gB^0 (4.16) 

If one deals, however, with non multiplicative operators, like e.g. differential operators, 
an additional commutator has to be added. Also note that unlike Tr , the symbol ( ) has 
dimension of . 

More explicit expressions for and (3^ are given in Table 4.2. Note that , 
depend only on D (mod 2). 
In particular for D — 4: 

/ 2^ 11 3 1 \ 

^ = ( 7T7^^(^o - ^nv'Oi -m'02 + 2mO^ - -O5 + ■ ■ ■ ) 

\ 24 9 2 m , , 

+ log (^] /-^m^Oo + Irn^Oi + rr^O^ - 2mO^ -o\ 

We can distinguish three contributions to the effective action, according to their m 
dependence, namely, 1) the contribution from oc^ and /5^1og// for N < D, which is a 
polynomial of degree D in m, 2) a logarithmic part from log m with N < D and 3) an 
inverse mass expansion from and N > D. The polynomial part is regularization depen- 
dent as it can be modified by adding suitable local polynomial counterterms to the action. 
These terms are such that the action depends on D and not on Dq and m separately and 
they vanish for odd dimensions. The logarithmic and inverse mass parts are regularization 
independent. In fact by applying the operator (d/dm)"^, n > D at both sides of Eq. (4.14), 
the l.h.s. becomes ( — l)"^Tr (D""^) which is ultraviolet finite and hence independent of the 
renormalization prescription. In the r.h.s. the polynomial part as well as the dependence 
on fj. disappears whereas the other terms remain, yielding a pure inverse mass expansion 
without logarithms. The dependence on log// was a trivial additive constant for the un- 
regularized action. This is no longer the case after renormalization, indicating that the 
action has developed a scale anomaly. The practical interest of the former arguments lies 
in the possibility of reconstructing the action by dealing with explicitly convergent and 
hence regularization independent objects, namely (— l)"Tr (D~"'), n > D, and integrating 
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back in the mass parameter. By properly fixing the arbitrary integration constants one 
might reproduce a given renormahzed action. This idea is aheady contained in Eq. (2.10) 
for the specific case of the ^-function. 
4.3 EfFective currents 

Finally, we can also obtain an inverse mass expansion for the current J introduced in 
Section 2. Rather than computing J through the variation of (On) in Eq. (4.14), we will 
use the closed expression (2.13). Recalling that (a;|D^|a;) has a simple pole at s = —1, we 
obtain the following expansion for the current 

J(^) = E [^N-D+P^N-D^og m^-^-i0^(x;Do) (4.18) 

where 



w 

(4.19) 

(^K = -scd+k{s - 1) = eK+i/?K+i 

Again we can check that = 0. Also comparing with the expansion for the action, 
one finds 

^ = (Sp7J^--«0„_i(x) (4.20) 

which is consistent with Eq. (4.12). A similar relation holds for the Seeley-DeWitt coeflft- 
cients [2]. It is instructive to try to obtain the same coeflftcients On{x) from a heat kernel 
approach. The main result is that in even dimensions, the Seeley-DeWitt coefficients of the 
second order bosonic differential operator — coincide with the inverse mass expansion 
coefficients of even order (92n(D) of the first order fermionic operator D. This result has 
interesting consequences and will be proven in Appendix B. 

5. Anomalies within the ^-function regular izat ion 

5.1 Chiral and scale anomalies 

As it is well known the effective action does not share all the symmetries of the classical 
action. In the (^-function regularization approach, this is because not all the symmetry 
transformations of the classical equation Y)(p — are also symmetries of the eigenvalue 
equation D0„ = Xn(t>n- Symmetries which are broken by a mass term classically, develop 
an anomaly at the quantum level. This is the case of chiral and scale transformations. 

As already mentioned under the combined set of vector and axial transformations 
only the axial ones present an anomaly, within a ^-function regularization. Vector gauge 
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symmetry remains unbroken. The axial anomaly is defined [15] as the variation of the 
effective action under an axial gauge transformation of D, that is 



D 



(5.1) 



where a{x) is an infinitesimal matrix valued function in fiavor space only. This transfor- 
mation is consistent with the general structure assumed for D, i.e. D = +Y, with Y a 
local function. Aa is the axial anomaly, which is ultraviolet finite and of dimension D in 
the external fields [15]. The best way of computing A a in the (^-function context is to use 
the same trick as in Eq. (2.12): 



Aa = -SaTt logD = -5a— Tr 

as 



— Tr ( s{m75,D} 



D 



s-l 



D 

s=0 



8=0 



= Tr (2ia75D^) := Tr (2ia75D°) 



(5.2) 



s=0 



Here we have used the trace cyclic property and also that (a;|D*|a;) is analytic at s = 0. 
Since formally D° would correspond to the identity operator, the previous result can be 
interpreted as a regularization of this operator. Note that the anomaly is independent of 
the scale Now using the expression for (a;|D'^|a;), n = 0, 1, 2, ... in Eq. (4.13) (for even 
D), we have 

Aa = {2iaj50D) (5.3) 

which is manifestly ultraviolet finite. With the expressions for O2 and O4 in Eq. (4.7) 
and after some algebra using only defining properties of the gamma matrices, we obtain 
simpler explicit forms for Aa in two and four dimensions 

(D = 2) Aa^ (2i«(x)75D2) 

{D = 4) 



^A = (2m(x)75( V + JdA^D)) 



(5.4) 



2 3 

These expressions are more easily arrived at by going back to Eq. (4.1) to compute 
Tr (2ia75D*)|s=o directly. For illustration purposes we quote in Appendix C more elabo- 
rated versions of Eq. (5.4) for a Dirac operator with the most general spinorial structure 
in two and four dimensions, after explicit evaluation of the Dirac traces. 

The scale invariance can be treated in a completely similar way in this regularization. 
The trace anomaly is the variation of the effective action under a scale transformation of 
D[l]. 

T> = i^+Y{x) — +e-'Y{e-'x), As = SsW (5.5) 

s 
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We can best compute the trace anomaly by rewriting the scale transformation as a 
homogeneous transformation for D 



(5.6) 



and using the (^-function regularization of the action 




s=0 



= e-^sTr ( D 
as \ 



D 



s-l 




(5.7) 



= eTr (D°) 



Likewise the axial case, the zeroth power of D appears in the final expression for the 
anomaly. Using Eq. (4.13) one finally obtains 



for the trace anomaly. 

From the last line of Eq. (5.7), the scale anomaly can alternatively be written as 



which is consistent with Eq. (4.14). Because the axial anomaly 6aW = Aa does not 
depend on /i, and using the cross variation condition [Sa, Ss] = 0, we conclude that As 
is chirally invariant. In addition it is also scale invariant. Moreover the scale anomaly 
vanishes for odd D because from Eq. (4.20) Sx{Od) = for an arbitrary variation. 

Another issue is that of the anomalous breaking of parity in odd dimensions and the 
related Chern-Simons action [31,32]. It will be considered in Section 7. 

5.2 Minimal form of the anomaly 

The presence of the anomaly indicates that the efi'ective action has terms which are 
not chirally invariant. As it is well known difi'erent regularizations in principle produce 
effective actions which differ in local polynomial terms. Given the fact that such different 
actions are related by counterterms, whose parameters are to be fixed anyway by some 
renormalization prescription, they are considered physically equivalent. These differences 
reflect in turn in the form of the anomaly. The different regularizations give anomalies 
differing by so called unessential terms, that is, terms which can be obtained as the variation 



As = eCo{OD) 




As = e/i 



dW 



(5.9) 
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of local polynomial actions. Bardeen [15] worked out the four dimensional case, including 
vector, axial-vector, scalar and pseudoscalar fields in the Dirac operator. He showed that 
the scalar and pseudoscalar fields did not contribute to the essential anomaly and moreover 
that the only essential terms were those of abnormal pseudoparity, i.e. containing the Levi- 
Civita pseudotensor and thus purely imaginary in Euclidean space. This Bardeen's or 
minimal anomaly were later shown to derive from the Wess-Zumino-Witten action [17,20], 
which of course is not a local polynomial. Since the work of Bardeen, it has been shown 
that nongauge fields [33,34], and internal gauge fields (i.e. transforming homogeneously 
under chiral transformations) [35] do not contribute to the essential anomaly. This is also 
suggested by the fact that if the formal variation of the action is regularized in a chirally 
covariant way, there is an obstruction to the integrability conditions which depends on the 
vector and axial gauge fields only [2,12]. 

Let us restate this result with our formalism in a way which is easily extended to 
higher dimensions and more general theories, i.e. relying only on algebraic transformation 
properties but not on the detailed coupling structure of the external fields. That is, let us 
show that in fact for a completely general Dirac operator in two and four dimensions in 
Euclidean fiat space one can write enough counterterms to bring the anomaly to Bardeen's 
form. To be precise, let D be of the form D = Dq + X where X is a local function, and 
Do and X transform independently under axial transformations 

Do — ^ e-^"T5Doe-^"^^ X — , g-^^^^Xe-^^Ts (5.10) 

A A 

Do contains the derivative part, and hence the corresponding fields transform inhomoge- 
neously. We want to show that all the contributions to the anomaly coming from X can 
be removed by counterterms. The construction is more easily presented by following the 
approach of Ref. [35] . Let D^ = Do + tX, then one can write the identity 

W(D)-W{no)^ [ dt5xW{JDt)^ [ dt /"d^a;tr(XJ(Di)) (5.11) 
Jo Jo J 

Here J(D) is the current introduced in Eq. (2.11). The axial anomalous contribution to 

the action containing X will be local polynomials if and only if the current has the form 

J(D) = Jc(D) + P(D) (5.12) 

where P is a local polynomial and Jc is a chiral covariant current, transforming as 
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so that 5aTt (XJc) = 0. Notice the opposite sign for the axial transformation as compared 
to the Dirac operator, Eq. (5.1). Indeed, the above decomposition of the current yields an 
analogous separation for the action, namely 

W{'D) = W{r>o) + dt J d^a;tr(XJc(Dt)) + ^ dt J d^a; tr(XP(Dt)) (5.14) 

The first term gives ^^(Dq), the second is chiral invariant and the last is the local poly- 
nomial counterterm. The observation that the current is of the form (5.12) was already 
made in Ref. [26] for chiral fermions in the presence of gauge fields. As already pointed 
out by Bardeen and Zumino, the fact that the total current admits such a decomposition 
is not obvious and requires a constructive proof for each case. 

Let us construct P explicitly for the two dimensional case. Clearly P must satisfy the 
conditions 

,,P.,,J.^, ,,P = (5.15) 

as a consequence of Eq. (5.12) and vector gauge invariance. Notice that since the anomaly 
involves the symbol ( ) and D is not a multiplicative operator (see the remark to Eq. 
(4.16)), cyclic property might not be applied in principle to compute SAa/^^- Nev- 
ertheless, the vector gauge invariance of the anomaly and the fact that ST) = X is a 
multiplicative operator, allows to do so, yielding 5xAa = (X{2m75,D}). One can check 
that the first relation is then satisfied by 

Po = --^(x|D|0) (5.16) 

Unfortunately this solution breaks vector gauge invariance. In order to construct P we 
should subtract a new polynomial Pi from Pq to reestablish vector gauge invariance. Also 
Pi must be axial covariant in order not to modify the already correct axial transformation 
of Po 

6gPi = 6gPo , Pi -j^ e*"TsPie*'*^s (517) 

The only object algebraically made out of D, transforming axially as Pi should, is 7^07^, 
which however is also vector gauge invariant and hence cannot match SgPo. Therefore we 
must resort to new objects or use more information on D to write enough counterterms. 
Let D be other Dirac operator with 

5gt) = 5gU=4, D -^e^^^^De^^^s (5.18) 
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If D has the standard hermiticity, (i.e. the Hamiltonian is Hermitian in Minkowski space), 
a solution is provided by D = — D'''. Now it is straightforward to obtain Pi by writing 
all the possible objects of second order with the correct axial transformation, constructed 
with D, D and 7^ and adjusting their coefficients to match ^^Pq. We finally find the 
solution 



In the previous construction Pq has been obtained by trial and error. In more compli- 
cated cases, the best way to proceed is to introduce a polynomial action in the Dirac oper- 
ator, Wo(D), from which Pq formally derives, since the number of possible terms decreases 
substantially. The fact that the anomaly cannot be subtracted by local and polynomial 
counterterms, prevents the existence of such an action in a literal sense. Nevertheless 
one can impose SxWo(D) = J d^a;tr(XPo(D)) and SaWo(D) = Aa modulo commutator 
terms, which would vanish if cyclic property were valid, i.e. if all the operators involved 
were multiplicative. That is what we mean by formal in this context. There are two key 
observations. First, that actions Wo(D) which are algebraically made out of D, uniquely 
determine the current Po(D), namely, by substituting D — > D+X in Wo(D), keeping terms 
with just one X and freely using the cyclic property to bring all the X say to the left. The 
relation is unique even if Po(D) is only the formal variation of VFo(D). And second, if the 
current Po(D) reproduces the variation of the axial anomaly, i.e. 5aPo(D) — 5.4.^/50 
(see Eq. (5.15)), the action must reproduce the anomaly 5^VFo(D) = Aa at the formal 
level. It is therefore advantageous to solve this latter equation and hence to derive Po(D). 

One can check that a solution in the two-dimensional case, is given by the action 



which of course is not gauge invariant. Wo(D) formally gives the axial anomaly, and 
its current is the same Po(D) found previously, Eq. (5.16). In summary, Pq is a local 
polynomial current with the same anomaly as J. Once Pq is available one can proceed as 
explained above to obtain P. 

The previous method can be applied to the four dimensional case. We find for Wq 




(x|(D-D)|0) 



(5.19) 




(5.20) 
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and P the following expressions 



^o(D) = -(-D^ + -D^^^D^^^ + -0^7^07^) 

^ = ~^ (4^^"^'^^^°' ' + D7^D%^ - D7^DD7^ 

+ 7mD^7m - 7mDDD7^ + D27^D7^ - DD7^D7^ 



(5.21) 



+ 7^D7^D2-7^D7^DD)|0) 



The terms without D are those coming from Pq. This polynomial generalizes that found in 
Ref. [26] . This completes the proof that all the fields in D transforming homogeneously un- 
der chiral rotations do not contribute to the essential anomaly in two and four dimensions. 
If one keeps only external vector and axial fields in Dq, ^^(Dq) is Bardeen's anomaly 
(up to some normal pseudoparity terms which can again be removed by counter terms). 
Explicit expressions for the counterterms after having worked out the Dirac algebra can 
be looked up at Appendix D. 

An interesting aspect of the previous results is that the counterterms needed to repro- 
duce the Bardeen form of the anomaly in ^-function regularization requires introducing, 
besides D, a new Dirac operator D transforming in the same way as the current under the 
chiral group. As already mentioned, an operator transforming in that way in Euclidean 
space is given by — D^. This agrees with similar findings in other regularization schemes, 



using ad hoc prescriptions like e.g. to separate the action into real and imaginary parts 
WCD) = W+CD) + W^-(D), by means of the formula 



and similar ones [2,10,36,37]. 

In contrast to the axial anomaly, the scale anomaly ^5 contains no unessential terms: 
by dimensional counting the possible local polynomials would be scale invariant or else 
would have to include external functions, thus introducing an anomaly in the Poincare 
symmetry. 

6. Gradient expansion 

6.1 General considerations 




(5.22) 



W^-(D) 



^{TV log(D2)-TVlog(Dt2)} 
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In Section 4 we considered an inverse mass expansion for the effective action of D. 
It was both an expansion in the number of external fields and the number of derivatives. 
Here we shall consider an expansion in the number of derivatives and the number of fields 
with Lorentz indices. That is, we take 

D = M + Do, M{x) = S(x) + i^5Pix) (6.1) 

where S{x) and P{x) are scalar and pseudoscalar fields and Dq includes as well as vector, 
axial vector, tensor fields, etc, and we expand in powers of Dq. This is a resummation 
of the inverse mass expansion, relaxing the restriction that S and P should be weak 
fields. Another standard resummation, complementary to this one, is the perturbative 
expansion which assumes weak but not necessarily smooth fields. The gradient expansion 
is a semiclassical expansion similar to that used in quantum mechanics and many body 
physics [38,39]. This means that its starting point approximates the spectrum of D by 
a continuum. Discretization effects are averaged and cannot be recovered in detail by 
resummation, hence the expansion is at most only asymptotic. The same is true of course 
for the inverse mass expansion. 

Expanding eq. (3.6) and using (3.7), we have 

oo 

^ ^ (6.2) 

Go{x,p; z) = i^-z + M)-^ 

= -{p^ + {z- Sf + P^+ i^,[S, P])-\i) +Z-S + ii^P) ^^'^^ 
Let us remark that ioT N > D the integrals are ultraviolet finite, yet we cannot proceed 
formally by simply taking (log(M+ ^ +Do)) and expanding in powers of Dq because 
this would require trace cyclicity and in fact gives wrong results (also we cannot expand 
formally Tr log(M + Dq) before using the Wigner transformation method because every 
term would diverge). However it can be done for the current or, equivalently, we can take 
s = — 1 above and perform the z integral first. Thus, we obtain for the current in A'^-th 
order, with N > D 

[JW]. = (-l)-'y^tr(.|(^D.) ^10), N>D (6.4) 
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The action can then be reconstructed with the Eq. (2.14). Because z no longer appears 
in Eq. (6.4), 3n>d enjoys all the classical symmetries of D; so long as the symmetry 
transformation does not mix different orders of the expansion, which is the case for the 
usual Poincare invariant internal symmetries. 

If Do contains only fields with an odd number of Lorentz indices, the terms with N 
odd vanish after performing the p integration in even dimension. 

6.2 Effective action in 1+1 dimensions 

The problem with gradient expansions is that it is in general difficult to work out 
explicitly the inversion of matrices implied in Go{x,p; z), which is necessary to perform 
the z and p integrals [9] . To be concrete and keep the computations simple, we will consider 
the 1+1 dimensional case with SU(2) flavor symmetry, that is 

D = z^ + f +475 + 5 + Z75P (6.5) 

where S{x) is a Hermitian c-number and P(x), V"^(x) and A^{x) are matrices in the 
fundamental representation of su(2). Because of the two dimensional identity 7^75 = 
i^ixvlv , we can simplify the amount of algebra without loss of generality by reabsorbing 
A^{x) in V^{x) with V^ = V^- ie^„A„. 

Let C be the matrices in the fundamental representation of sl(2,C), i.e. of the form 
ttiTi where Tj are the Pauli matrices and are c-numbers. The following properties of C 
will be widely used later: 

for A,B e£, tvA = 0, [A, B] e £, {A, B} = c-number (6.6) 

In particular Go{x,p; z) can be computed to obtain 

where we have defined 

a = z-S A = p2 + (72 + P^ iV^ = id^ + % (6.8) 

and o", and A are c-numbers. Odd orders vanish. Let us compute explicitly the first 
two terms. A?" = 0, 2. One finds integrals of the following form 

Uk,N,D,s;S,P) = j ^ j^^^ (6.9) 
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Because the integral is a meromorphic function of D, we can perform first the integral in p 
as in dimensional regularization, and then the z integral gives the hypergeometric function 
2F1 [40, p. 555]: 



(47r)^/2 r(f)r(iv)r(n-s) V ' ' ^+ 

+ e-»'^'^zi-^"+^ 2^1 (n, 2n - s - 1; n - s; ^) 

(6.10) 

with n = N — |D — A;, z± = S ± i\P\ and |P| = {P'^)^. For even D it is convenient to 
write the same result as 

1 T{N-n) / ^ 

2 - (47r)i^/2 r(f )r(iv) [d^J {z - ^' > 

with the prescription 

In this form eq. (6.11) holds too if z^ is replaced by any function analytic in C— R~ and 
sufficiently convergent at 00. After taking the trace in Dirac space we have 

And hence for the action at zeroth order 

.2 A 



Wo = i^M' (^log(-^) - 3j + 2P'^ (6.14) 

( ) was defined in Eq. (4.15) and in these formulas tr acts on flavor space only, and 
— S'^ + P'^. A similar calculation in four dimensions gives 

Wo = - (m^ (log(^) - f ) + \m^P' + ^P'^ (6.15) 

The polynomial chiral breaking terms yield the spurious axial anomaly discussed in the 
previous section. 

Let us compute the second order term. After taking the trace in Dirac space, we can 
always proceed by carrying the covariant derivatives iV^ to the right, for instance 

zP^A-i = -A-\-2a[iT)^, S] + {P, [iT)^,P]}) + A-hV^ (6.16) 
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so that the various A and a are collected to the left and the integral I2 applies. This is 
not the most efficient strategy in this case but it is systematic and allows for algebraic 
manipulator implementation in more complicated cases. Due to vector covariance, one 
finds that ib^ appears only covariantly, that is in the form [iP^, X]. The vanishing of the 
non covariant terms follows after the p integration and using the C identities, eq. (6.6). 
The result can be brought to a more symmetric form by using integration by parts of the 
covariant derivatives and trace cyclicity: 



^ ^ 12 2 M2 

P - lliS"! P ~~ 

+ arctan(-)ie^^F^^ + -^{^j^ - ^BxciQxv{-))e^^PD^PD^P 



(6.17) 



Here D^X = [V^,X], iF^^ = [W^,,W^]. 

6.3 Analysis of the effective action. Wess-Zumino-Witten action 

In order to obtain a separation 14^2 — + in two terms with well defined 
hermiticity and pseudoparity, we have to reexpand iV^ as -D^X = D^X — e^,^[Ai,, X]. The 
real part, W2~, comes from terms without explicit e^j^ and even number of axial fields or else 
with explicit e^,^ and odd number of axial fields, whereas W2 is the other way around. In 
principle, the effective action up to second order in the gradient expansion should saturate 
the anomaly equations. This can be trivially checked for the scale anomaly. Although 
the resulting expressions are manifestly vector gauge invariant, the corresponding check 
for the axial anomaly requires more work. A convenient method to do so is to explicitly 
establish that the anomalous terms of the computed effective action coincide with the 
Wess-Zumino-Witten action plus additional polynomial counterterms. 

To analyze the chiral transformation properties of W2, let M{x), U{x) and (j){x) be 
defined by 

S + iP = MU, U = 6'^^ = cos + i sin (j) (6.18) 

hence U G SU(2), sin(^ G su(2) and M, cos(/> are c-numbers. (f){x) is a multivalued 
function of S and P, well defined everywhere by assuming explicitly that M{x) > 0, or 
equivalently that M is non singular, as required by Seeley's formula (2.7). 

Under the previous assumptions, the terms in W2 with explicit e^i, turn out not to 
depend on M{x), as it is readily shown. That is, this part of W2 can be written as 

^1^- ^ / 2iF^ + 1 ^^^(2*^) - 2^ sin0(^sin(^)A (6.19) 
\ 4 sin / 
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i^ = (f) + rnij^, neZ (6.20) 
\4>\ 

(f is such that tan(^ = tan0 and takes into account exphcitly the multivaluation. A 
standard exterior algebra notation has been used [26] , 

F = ^F^^dx^dx^ , DX = D^Xdx^ (6.21) 

and dXfj, are anticommuting variables. For more general symmetry groups, M{x), defined 
by the factorization S + iP — MU with M Hermitian and U unitary, will not be a c- 
number and will not cancel in W~ . Nevertheless, the symmetry under local rescaling 
M(a;) A(a;)M(a;), where X{x) is a c-number local function, still holds. This general fact 
will be established below when dealing with the fermionic current. 

In order to work out the expressions it is convenient to use the following identities, 
particularly useful for computing derivatives of the fields: 



foraeA SU{a) = ^{aJa}, Sf_ia) = ^-^{a,5a} + ^[a,Sa] (6.22) 



where f± are arbitrary even and odd functions respectively, and /4 their derivatives. As 
a consequence 



1 1 

iuHu = -- — -5 sin H — (tan + [sin 0, 5 sin 0] 

cos (p 2 

1 1 

iU 5U^ = H -5 sin (tan (j) — i) [sin 0, 5 sin 0] 

cos 2 



(6.23) 



After some algebra, the real part of W2 can be written as 
where we have defined 

= uHd^u - + u^v/^u , F/'^ ^v^±A^ (6.25) 

7?.^ transforms covariantly under chiral gauge transformations 

^r,l{^) e SU(2), U ^ QlUQ^j, , y^^'^ ^ ^^i^,L(^a^ + "t^^''''')^^J^,L (6.26) 

Vector gauge transformations correspond to the diagonal subgroup. M is chiral invariant 
and the only chiral breaking term in W2 is which is a polynomial. 
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Similarly, the imaginary part can be brought to the following form 

= rwzw(t^) + TaiU, V^, V^) - Toil, V^, V^) 
TwzwiU) = I HUhdUf (6.27) 

where dX = d^Xdz^. Fwzw is the correctly normalized Wess-Zumino-Witten action 
in two dimensions. Tg{U,V^,V^) is obtained by chirally gauging Fwzw [20]. The last 
term is known as Bardeen's subtraction, a counterterm which reestablishes vector gauge 
invariance. There / stands for the identity of SU(2). 

The way Fwzw appears deserves some comment. The field configuration U{x) defines 
a two dimensional manifold M2 without boundary inside SU(2), and D3 is a three dimen- 
sional manifold such that dD^ — M2. Up to a quantized multivaluation [20], Fwzw does 
not depend on the choice of D3 because ^3 = j^tr{UHdU)^ is closed. Hence, there is a 
2- form Q2 such that locally ivs = dfl2- However ^2 cannot be regular everywhere, because 
(jUs is not exact (in fact it is just the volume element in SU(2)). Choosing D3 to avoid the 
Dirac singularity one has 

Fwzw= / a;3= / ^2 (6.28) 

To arrive to eq. (6.27), what one can show is that a valid choice for VI2 is 

^2 = -^tr / ^in(20)-2,. ^.^ \ 
167r \ sin / 

which is just 1^2'"^ (6.19) for = = 0, and arbitrary n. The singularity is at 

U = (—1)'^'^^. In fact the solution of ujz = dfl2 is unique imposing global vector invariance 
and regularity at = 0. Global chiral invariance is not manifest for any choice of Cl2- Also 
one can check that by varying n, J Q2 changes by an integer multiple of 27rz. Another 
interesting point is that vector gauge invariance would be achieved by minimal coupling 
zd — > zd + F in il2- However such an action would not be single valued modulo 27rz. The 
new vector gauge invariant term {2iFip) in W2'"' reestablishes the one-valuedness of the 
action. 

can be written in other interesting form by applying Stokes' theorem to Fg, 

namely 

W2=^f tiiln^ - iUTr - 2%ATr + \a^) - p.c. (6.30) 
47r Jdz 6 6 
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where TZ = UUdU — Vr + W^VlU, Tr = dVji — iV^ are chiral covariant, and p.c. means 
to subtract the parity conjugate terms, i.e. the same terms exchanging U ^ U''' , Vr Vl 
and A <-> —A. This form is manifestly vector gauge invariant and the axial anomalous 
terms (those with A^) are independent of U. These terms are polynomial, yet they cannot 
be removed by counterterms because they do not form a closed 3-form by themselves. 

We would like to emphasize that ours is an ab initio calculation of the action. Once 
the ^-function prescription is adopted, there is no more freedom nor ambiguity in the 
calculation. This is in contrast to derivations of F-wzw by integration of the chiral anomaly, 
which have to assume that S + iP lies on the chiral circle, i.e. M{x) constant [41-44]. 

6.4 The fermionic number current 

As another illustration of the Wigner transformation technique, let us consider the 
fermionic number current —ip^nip in four dimensions with arbitrary internal GL(n,C) sym- 
metry. It is obtained from 

= - (x), S^W = J d''xuj^ix)j^ix) (6.31) 

where u!fj,{x) is a c- number field. As an example, for D = z ^ -|- ^ (x) + m, a, simple 
calculation from Eqs. (4.7) and (4.17) gives 

W = -l{F',,)log(^^^+0{m') (6.32) 

which gives a contribution to jn{x). Here we will consider the anomalous part of the 
current, jj^{x) in the presence of scalar and pseudo scalar gl(n,C)-matrix valued fields 

D = z ^ +M(a;) , M.{x) = S{x) + i^5P{x) (6.33) 

Because the pseudoparity odd (p.o.) part of the current is ultraviolet finite we will use 
directly the expression 

6^W- =Tj: (6.34) 

On the other hand, the lowest order contribution appears at fourth order in a gradient 
expansion of W, due to the presence of the Levi-Civita pseudotensor. Thus to lowest order, 
the anomalous current is obtained from 
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It is convenient to introduce the definitions 



MR^L = S±iP, F 
Glr = ip" + MlMr) 



-1 



±75) 

Grl = ip" + MrMl) 



-1 



(6.36) 



so that 



-1 



+ M}M)-^{i) -M)) = 
PrGlrMlPr + PlGrlMrPl - PrGlr ^Pl- PlGrl ^ Pr 



(6.37) 



Substituting in J^^VF , and keeping only terms with 75, one obtains 




(6.38) 



The parity conjugate (p.c.) is obtained by exchanging the labels R and L everywhere, 
tr no longer includes Dirac trace, and an exterior algebra notation has been used with 
Lo = u)fj,dx^ and d = dx^d^. d derivates to the right until it finds another d or gives zero if 
it reaches the right end. The integral over p is convergent, but it cannot be done in closed 
form for arbitrary Mr^l fields. 

An interesting property of Su,W~ , which can already be derived without explicitly per- 
forming the momentum integral, is its invariance under the local rescaling S{x) — > X{x)S{x) 
and P{x) X{x)P{x), with X{x) an arbitrary local c- number function. Due to the close 
relation between the fermionic number current and the imaginary part of the effective 
action [45,46], to be addressed below, this property holds for W~ at lowest nonvanish- 
ing order as well. As a consequence, in the particular case of Mr{x) = M{x)U{x), 
M{x) being a Hermitian c-number, not necessarily constant, the expression simplifies: 
Glr = Grl = (p^ + M^)~^ is also a c-number, all the terms containing dM{x) are readily 
shown to cancel, the terms with vanish and the integral over p is immediate, with the 
following result 





(6.39) 



This gives the correctly normalized Goldstone-Wilczek current [19] 




(6.40) 
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Comparing Eqs. (6.27) and (6.40) we find that the 3-form us is the integrand of the 
action in two dimensions and also it is the fermionic number density in four dimensions, 
in both cases at fowest order and in the pseudoparity odd sector. 

6.5 Effective action from the fermionic number 

The relation pointed out in the last paragraph of the previous subsection is much more 
general as noted by several authors [45,46]. We will show subsequently that our formalism 
can handle the proof of this general relation quite naturally. The general statement is as 
follows. Let D be a Dirac operator in D dimensions and W~ the pseudoparity odd part 
of its action. Consider the extension of D to a family of Dirac operators D(t)), v G [0, T] 
which interpolates between D at v = T and Dq at v = 0, where Dq will be characterized 
below. Similarly consider the extension of the space-time with two more dimensions {u, v), 
u playing the role of Euclidean time, and of D(f ) to a -D + 2 dimensional Dirac operator 
D': 

D' = ^uidu + ividy + D 

(6.41) 

7^ = X CTi , 7„ = Z/ X (72 , 7^ = 7m X CTs , 75 = 75 X (73 

The D + 2 matrices 7^ satisfy the Dirac algebra {7a) 76} = —'^^abi a, 6 = tt, 0, . . . , D — 1, 
and 75 = i'ju'lv^b for even D. The operator D is defined by making the substitutions 
7/i,5 ~^ 7m, 5 ill D(v). Thus it has the form 

where D''" appears because the hermiticity of the fields in D is such that the Hamiltonian 
luilvidv+fy) of the D+2 dimensional extended Dirac operator, is Hermitian in Minkowski 
metric, and hence D' has the standard hermiticity in Euclidean space. 

The theorem establishing the relation between the fermion number and the effective 
action can be stated in the following manner. Let B be the D + l-dimensional spatial region 
defined by {v,x^), v G [0, T] and G R"^, and F the fermion number enclosed in B in 
the D + 2 dimensional system. Then for a suitable choice of the interpolating path D(z;), 
to be considered below, 27rF = —W~ in Euclidean metric, W~ being the L>-dimensional 
pseudoparity odd effective action. In fact, one can state this relation in a more general 
form, namely 

J^d''xdvtr{u,v,x\^J'{t)')\u,v,x) = ^ (IV /(D) -Trf{D^)) (6.43) 
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for a sufficiently convergent function f{z) analytic in C— R~, and f'{z) its derivative. 
The relation between F and W~ is obtained with f{z) = log 2; defined as the analytical 
extrapolation of d2;*/ds in s at s = 0. 

Let us consider under which conditions this relation will be valid. The contribution 
from the value v = to the r.h.s. of Eq. (6.43) will vanish if Dq is related to by a 
similarity transformation. This can always be achieved for even D by taking Dq without 
pseudotensor fields, i.e. without explicit 75, since in this case Dq = 75D075. Observe that 
D' is time independent and does not contain fields with Lorentz index u, hence u in the 
l.h.s. must come from the Levi-Civita pseudotensor and it is automatically pseudoparity 
odd as the r.h.s. Regarding the choice of the interpolating path D(f ), note that the l.h.s. 
can be obtained from 

= - , 5^Tr /(D') = J d^x du dv tUaJa (6.44) 

Where c/= '^u^u + Jv^v + T^i^^^i and fx = 1, . . . , _D. Because of vector gauge invariance, 
Ua can only appear in the form daUJb — dbUa, hence locally = dbJab with Jab anti- 
symmetric. Note that the l.h.s. of Eq. (6.43) can also be written as — Jg d^xdvJu = 
— J^d^xdv{dyJuv + d^Ju^)- The terms containing J^^ cancel because we are assuming 
D boundary less in R'^. The non vanishing contribution comes from J^v at the boundary, 
i.e. V — and v — T. However, for even dimensions the v = contribution vanishes 
since Juv must contain a, D + 2-dimensional Levi-Civita pseudotensor, which cannot occur 
if Do is chosen as explained above. If we take an adiabatic path, that is, a smooth path 
with dy'D{v) = 0{1/T) and T — > 00, Juv\v=t will depend only on D and f{z) and not 
on the particular adiabatic path chosen. For arbitrary paths J^v ait v = T would also 
depend on 5"D|„=t, n > 0, hence the result would be path dependent. Let us remark 
that at lowest order Jab exists only locally if f{z) is a multivalued function such as log 2;. 
In this case different homotopy classes of adiabatic paths pick different branches of f{z). 
An example is again provided by UJ3 which is closed and hence locally exact. Indeed one of 
the interesting aspects of the relation 27rF = —W~ is that F is computed as the integral 
of a single valued expression. 

To proof Eq. (6.43) let us rewrite the l.h.s. using Seeley's representation and the 
Wigner transformation trick in the space spanned by \u,v) 

-Tr [ dvf^-^ [ —f'(z){u,v\% . 10) (6.45) 

Jo J (27r)2 Jr2m' ^ ' ' %{idu + cuj + ^,{id, + a;,) + D - ^ ' ^ ^ ' 
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|0) is the state with zero momentum in the space |tt, v), i.e. {u,v\0) = 1, and and dy 
act to the right until annihilating |0). Tr includes trace in the space \x^) as well as spinor 
and internal degrees of freedom in D + 2 dimensions. The ultraviolet limit corresponds 
to integrate over uj which is finite after the z integration. To proceed let us introduce the 
definitions 

r = i%% = ixa3, r = rt , = / 

(6.46) 

iH = %i-D-z), {H,T} = 
Manipulating the denominator of Eq. (6.45) as in (6.37) the formula can be written as 



[idu + cou) + + '^v) + -n - iT{dyH) 

(6.47) 

We can set du = everywhere and also dy\0) = 0. Furthermore in an expansion in powers 
of dy, one can check that only odd powers survive after doing an angular average over iv 
and taking trace using that H is off-diagonal. In particular for an adiabatic interpolating 
path, only the term linear in dy survives. Then the expression takes the form 

Now we have the trace of an operator in the space |a;^) only. Because f{z) is convergent 
the trace exists and the cyclic property can be freely used in the two terms in Eq. (6.48). 
The first term can be put in the form dy{Tg{H'^,(jj'^)), and it vanishes after taking the 
trace because H'^ is block diagonal with blocks (D''" — z) (D — z) and (D — z) (D^ — z) which 
have the same spectrum. In the second term we can integrate uj, thus 



_J_Tr [ ^f'(z)TdyI)J— = ^ r dvdyTr{Tf{i})) 

Jo Ti-z 47rJo (6.49) 

'^(r/(D))i:=o 



471 



which coincides with the r.h.s of Eq. (6.43). This completes the proof. Let us remark that 
the staircase relation (6.43) is valid beyond a gradient expansion (in D dimensions), and in 
case this expansion applies the leading terms provide a relationship between a topological 
action and a topological fermion current. Furthermore Eq. (6.43) does not rely on the 
chiral circle constraint and holds for any space time dimension. 
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A rather amusing illustration of the relation comes from the zero dimensional problem. 
In D = the most general Dirac operator is D = S + i^^P, S, P being Hermitian matrices 
in flavor space and 75 a 1 x 1 matrix which equals —1 in our convention (see Eq. (2.2)). 
Certainly, the Grassmann integral (2.5) can be carried out and W = — trlog(D///) exactly, 
without regularization. Remarkably, the effective action is ultraviolet finite yet it possesses 
both scale anomaly and (essential) axial anomaly (even if there are no gauge fields). This 
is also in agreement with the general formulas (5.3) and (5.8). 

In the Abelian case D = Me^'^'^'^, M, (p GR, the action is simply W = — log(M/|U) 
On the other hand, the pseudoparity odd singlet current in two dimensions (0 + 2) is 

= -^^f^'^dM + ^9lcp +...) (6.50) 

With fiji/ja = u, V. In the notation used in the proof of the theorem stated above, this 
formula corresponds to identify Juv = + ^v^/i^^^) + •"■); since does not 

depend on u. As expected the fermionic current at lowest order is a closed 1-form whereas 
the higher orders are total derivatives of v. j~{x) checks the general Eq. (6.43), the 
higher orders vanishing only in the adiabatic limit. Note that in the non adiabatic case 
the fermion number depends on the v-derivatives of at v = and v = T, and hence does 
not agree with the zero dimensional pseudoparity odd effective action. 

7. Odd dimensions 

In this section we will consider the odd dimensional case in greater detail. The main 
issue here is the quantum realization of parity which is a symmetry at the classical level [31] . 
For odd D and up to a similarity transformation, parity corresponds to the transformation 
D = +Y{x) — Y(— x). It turns A„ into — in the eigenvalue equation (2.4). As 

a consequence we can expect that divergent terms in the action can break this symmetry. 

In an inverse mass expansion of the action we must look for invariance under m — > —m 
and Do — > — Dq, however the formula (4.14) is not appropriate because it was derived for 
m > 0. The cut F = R~ in the integral Ii breaks parity explicitly, and indeed for negative 
m the free spectrum m ± i\p\ would intersect F. In fact, there is no ray of minimal growth 
for both positive and negative masses as required from Seeley's representation, Eq. (2.7). 
To deal with this problem one can use some infrared regulator such as including a factor 
d{\p\ — mo) in the p integral and letting mo — > O"*". However, a simpler procedure is to use 
dimensional regularization instead of ^-function, which also maintains gauge invariance. 
As always, both regularizations will coincide up to counterterms. The procedure in this 
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case is to compute the current J = — (a;|D~^|a;) in a formal inverse mass expansion for a 
generic dimension and afterwards to reconstruct the action. Of course we do this using 
the Wigner transformation. One easily obtains the result 

J(^) = TX^^™)"" E m^-''-'Oj,{x; Do) (7.1) 

^ ^ N=0 

where D = Dq + m, 0{x) are the same coefficients as in Section 4, and 



TTl 

e(m) = — r, m e R 
m 



a 



1 , N 



(7.2) 



2 "2 

The terms of J in (7.1) which are ultraviolet finite coincide with (4.18) (including the factor 
e{m)^ which comes from integral Ii by making the momentum integral first). On the other 
hand, the terms N < D oi the current depend on the regularization. The coefficients a}v<£) 
have poles at even integer values of D. Such poles give rise to the logarithmic terms in 
(4.18). For odd D, the terms N < D in (7.1) give a finite contribution proportional to e(m), 
whereas in ^-function they vanish for positive m. A more careful treatment would show 
that in fact there is a step function factor 9{—m) in (4.18), so that both regularizations 
differ only by a polynomial of m, since m^9{—m) = m'^l — ^m^e{m). The sign function 
e(m) does not necessarily vanish at m = 0; its precise value is ambiguous and further 
information is required. 

The current in (7.1) is manifestly gauge covariant and also parity covariant (for 
odd D) thanks to the factor e(m), since both symmetries are respected by dimensional 
regularization and the current turns out to be finite. In contrast, to compute the ac- 
tion directly in dimensional regularization one would need a Seeley-like representation 
Tr log(D) = — (27rz)~^Tr j-^dzlog{z)(D — z)~^, thus manifestly breaking parity invari- 
ance. This situation resembles the method of Leutwyler [12] and Ball [2], based on the 
integrability obstruction of the covariant regularization of the formal variation of the ac- 
tion. We must then consider whether J is a consistent current and independent of the 
choice of the splitting of D into Dq and m, that is, whether dmJ — 0. To answer these 
questions let us consider separately the ultraviolet finite and divergent parts, 3 — 3 f + 3(i, 
respectively. The finite part coincides with the ^-function current and hence it is consis- 
tent, indeed it is the variation of the action (4.14) (after including the e(m) factor). Also 
it is invariant under Sm by the same arguments as in Section 4, with the prescription of 
defining e{m)m^~^'^^ as a derivative of the distribution e(m) log |m| in the m variable. 
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Consider now the divergent part of the current, corresponding to N < D, and odd 
D. From (4.12) it is immediate to check that Sm^d = ii m ^ 0, that is, except for Dirac 
delta terms coming from Sm^{'m). Hence we can resum all these terms by pulling out the 
factor e(m) and setting Dq — > D and m ^ in the remainder, thus 

J'^ = -^4^^--i(-;D). (7.3) 

Because (Od) vanishes for odd D, it is not obvious whether is consistent. We can check 
that up to the e(m) factor, to be discussed below, the current appears to be consistent 
for D = 1 and 3. Indeed, this can be accomplished by writing the possible terms in the 
action and adjusting their coefficients. In the case D > 1, one can see that this cannot be 
done with just terms of the form (P(D)), P being a polynomial, and one must resort to 
polynomials on D = +Y and T) = — Y. For given J the solution is unique up to an 
absolute constant. We find 



(D = l) Wd = -'-^ J dxtT{Y) 

(D = 3) Wd = I d'xtT (1y^ + ^Y7^Y7^Y + ^-Yi^ Y) 



(7.4) 



In the D = 3 case 707172 = cr = ±1 distinguishes the two inequivalent irreducible repre- 
sentation of 7^, and Y =4 +M, M being Hermitian flavor matrices. Working out the 
Dirac space algebra, the divergent parts of the current and the action can be written more 
explicitly as 

Jd = ^(2M2 + lae^^aiFf^^la) 

Stt 2 . 

e(m) f 4 2 

Wd = / d^xti (-M^ - iae^^a{A^d^.Aa - -iA^A^A^)) 

Here tr no longer acts on Dirac space. Up to the factor e(m), the pseudoparity odd term 
of the action, , is the correctly normalized Chern-Simons action [31] and it is gauge 
invariant. 

Finally let us consider the issue of whether is invariant under 5^, i.e. whether it 
is independent on the choice of the expansion point m. The invariance is broken by the 
factor e(m) in (7.3), and Sm^d = —{4'n)~^^~^^^^OD-i8{'rn)5m does not vanish as a distri- 
bution. Still one can consider adding counterterms to cancel the variation. The rationale 
is as follows. The ultraviolet divergent integrals become well defined by applying d/dm 
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a sufficient number of times, hence the ambiguity in the integral must be a polynomial 
in m. By the same token, expressions which diverge as m — > become well defined after 
multiplying by m a sufficient number of times, hence the infrared ambiguities, i.e. the 
admissible counterterms, consist of the distribution 5(m) and its derivatives. Unlike the 
ultraviolet case, however, there are infrared divergences to all orders in the inverse mass 
expansion. 

For instance in the case D = 1, if we change = — |e(m) to = — |(x|e(D) |0), and 
expand in powers of Dq, the difference are just infrared counterterms and now 6mJd = 
explicitly. Unfortunately this current is no longer gauge invariant nor consistent. So wc 
must consider simultaneously the three conditions of Sm invariance, gauge invariance and 
consistency of the current after introducing infrared counterterms. Considering again the 
D = 1 case one can see that these conditions are not compatible. Indeed let 70 = icr, 
Y = M + 7o^o> M,Aq gR. 5m invariance implies that 3 a must depend only on the 
complex number Y, consistency requires that 3d must be an analytic function of Y and 
finally gauge invariance requires J to be independent of ImY, therefore must be a 
constant, which contradicts = — |e(m) for Y = m. Also, this solution would break 
parity, which requires the current to be an odd function of Y. 

In conclusion, if we insist that the current and the action must be gauge invariant 
we must choose between two possibilities. First, take m and Dq as independent variables. 
In this case the actions (7.4) are simultaneously gauge and parity invariant. And second, 
that the action depends only on D. In this case we obtain two different actions by taking 
m positive or negative both of which are gauge invariant but have a parity anomaly. This 
is what happens also in the massless case [31]. The choice between the two possibilities 
depends on the physical problem at hand, for instance whether the problem admits a 
natural definition of m or not. Note that the ultraviolet finite part of the action was 
independent of m. 
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Appendix A. Seeley-DeWitt coefficients for Dirac operators to all orders 

In the case of even order, the Seeley-DeWitt coefficients for D admit a simple form. 
The general pattern can be illustrated by inspection of the first coefficients up to Oq(D) 



Co(D) = 1 
02(D) = + A^ 

C»4(D) = -D^ + ^(D^A^ + AD^A + A^D^) + -A"" 
2 3 6 

Oe{JD) = -D^ + — (D^A^ + D^AD^A + D^A^D^ + AD^A + AD^AD^D^A^D^) 

6 12 

+ — (D^A^ + AD^A^ + A^D^A^ + A^D^A + A^D^) + — A^ 
30 ^ ^ 90 

(A.l) 

Here, each monomial of degree 2m in A, stands for the (2m — 1)!! terms obtained by all 

possible contractions in the Lorentz indices involving the 2m A^ vectors. (Compare with 

Eq. (4.7) for O2 and O4). Note that does not commute with the symbol A and also 

that the cyclic property cannot be applied. 

One realizes that these coefficients depend only on A^ and D^. This is in fact a 

general rule which can be inferred from a direct comparison of inverse mass and Heat 

f n -\- 2m \ 

Kernel expansions (see Eq. (B.7)). Moreover, we note that the I monomials of 



n 

degree n in and 2m in A, with 2n + 2m = N , appear in all combinations and with the 
same coefficient ajv.n; namely 

In summary, the final formula for Cjv(D) reads 

N/2 

O^r(D) = J2 aAr,n[D2"A^-2-] (A.3) 

n=0 

Where [D^^A^-^^] stands for the sum of the terms obtained by all possible 

monomials of degree n in and N — 2n in A. As expected, the coefficients do not 
depend explicitly on the space time dimension D. In fact, if one assumes the validity of 
Eq. (A.3), the coefficients can be completely determined by taking particular cases. For 
instance, the relative weight of the coefficients is determined by gauge invariance, and the 
global normalization can be fixed by either going to the D = case or by taking D to be 
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a c-number. Finally, we mention that the odd order coefficients can be derived by direct 
use of the formula SmOpf = — e^v-DC^jv-i'^"^. 

The previous result can straightforwardly be applied to compute the chiral anomaly 
in D dimensions for a purely vector external field. Defining the field strength tensor 

iFfj^t, = [iD^.iD^], one finds 

where ( ) was defined in Eq. (4.15) and tr includes flavor space only. 



Appendix B. Inverse mass and Heat Kernel expansions for Dirac operators 
in even dimensions 

In even dimensions a straightforward relation between the well-known Heat Kernel 

expansion and the inverse mass expansion of Section 4 can be established. Our result is 
a reminiscent of the formal determinantal relation Det (D)Det (— D) = Det(— D^). Ac- 
tually, what is found is that this relation holds for the ^-function regularization of the 
determinants at each order in an inverse mass expansion, where the l.h.s. stands for the 
Wigner transformation result and the r.h.s. for the Heat Kernel expansion. 
Let A be a positive definite second order differential operator of the form 

A = -Vl + Y{x), iV^^id^ + B^{x) (B.l) 

The heat kernel {x\e~'^^\x) is ultraviolet finite for r > and it admits an asymptotic 
expansion [29] around r = 

oo 

(a;|e-^^|a;) = (Attt)-^/^ ^"«n(a^) (B.2) 

n=0 

where the Scclcy-DeWitt coefficients an{x) arc covariant local polynomials in iV^ and Y 
of degree p, q respectively with 2p + q = n. This coefficients can be written in a dimension 
independent way [2]. If we proceed formally, relying on the good behavior of the required 
analytical continuations, we find 

^ ' ' ^ r(-f) 7o \ I I / (B_3) 

A = -D^ + m? = -Do - 2mDo 
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Now, by using the heat kernel expansion, the r integral can be carried out order by order 
(for < and then analytically continued to > 0). In this way we obtain another 
representation of the ^-function 

(xlD'lx) = 5^^(-l)*-"!^-^-L^a„W (B.4) 
This is to be compared to the expansion obtained from the Wigner transformation approach 

D+8-N _|_ 1^ 

The comparison however is not immediate because we must first of all reexpand an{x) in 
powers of m, or equivalently in powers of and X, where Dq = +X. For a given 
order A^, On{x) generically gets contributions from all the terms an{x) with n > ^N. 
This is because an{x) is made out of and Y{x) which are of zeroth order in iV^^ and 
X, namely 

zP^ = iVf^ - ^{7^, X + m} 

1 1 1 ^^-6) 

y = - (X + m)2 - -[7^, \iV^, X]] - -{7^, X + m}2 - iV,] 

Nevertheless, whenever Ojv(a^) requires only a finite number of heat kernel terms we can 
check that both expressions coincide. In particular choosing D even, s = m = 0, we find 
the following identity {A7r)^/'^{x\'D^\x) = Od{x;'D) = aD/2{x;A = -D^), which holds 
regardless of the way the coefficients are written. If both coefficients are expressed in a 
dimension independent form, one has the further identity, for even dimensions 

02n{x;'D) = an{x;A = --D^), n = 0,1,2,... (B.7) 

This can be better checked for lowers orders rewriting Eq. (B.6) (for m = 0) as 

iV^ = -A^, y = -D^ - Al (B.8) 

Using well known expressions for ao, ai, 02 [2], we reproduce Oq, O2 and O4. An interesting 
consequence of Od{x;D) = 0^/2(0;; A = — D^) is that, although an anomaly calculation 
using as regulator might be questionable in principle for non normal D or non positive 
— D^, such a procedure turns out to be justified a posteriori. 



40 



Appendix C. Explicit form of the £Lxial anomaly in four dimensions 



The most general spinor structure for the Dirac operator in four Euchdean dimensions 
is given by 

T> = 11^+8 + i'y5P+ + ^ia^.T^. - ^ct^.TsT;, (C.l) 

with iVfj, = idfj^ + V^, cx^,, = lbiJ-,liy] and T*,, = \e^upaTpa. If the Hamiltonian 7°D 
is Hermitian in Minkowski space, and assuming that the fields transform as usual under 
the Wick rotation, the terms 5, 4 75 and \i(T^^Tp_i, are Hermitian whereas z^, ^75^ and 
icr^jy75T^j^ are antihermitian. This allows to distinguish between the two tensor terms. 
However, as our expressions will not depend on the hermiticity properties of D we can 
reabsorb T'^^ into T^^, i.e., T^^ - iT'*^ T^^. 

In the following expressions, ( ) stands for J d^x tr ^^^^(a;) and tr refers to flavor 
only. The axial anomaly is obtained adding up all the terms. 

A[V^A''] = -ie^^pA ^{\iV^,iV,],ApA^} + '^A^\iV,,iVp]A, - ^\iV^, A,]\iVp, A,] ) 
A[V^A] = i{ \iV^, \iV^, \iV,, A,]]] - \iV^, \iV,, \iV^, A,]]] + [iP^, \iV^, A,]]] ) 

AiVA"] = i{ [iV^, A,]A,A^ + A^A,[iV^, A,] + ^A,A,\iV^, A^] 

2 2 2 

+ -A^\iV^, A^]A, + -A^[iV^, A^\A^ + A^\A^, ) 

A\SP^\ = ( P{S, P}P + {S, ) 
A[SPV^] = { {[iV^, P], \iV^, 5]} + \{\iV^, [zP^, 5], P} + {[zP^, [zP^, P]],S} ) 
A[S^VA] = i{ {S\ \iV^, A^W + SA^\iV^, S] + \iV^, S]A^S 

+ ^[zD^, S]SA^ + \a^S\iV^, S] + ^S\iV^, S]A^ + ^^^[^^m: S]S ) 
A[P^VA] = i{ [[A^, P], [iV^, P]] + {[iV^, A^],P^} ) 
A[SPA^] = { [A^A^, [S, P]] + {{A^, S}, {A^, P}} ) 
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A[P^T^] = \i{ t^uP^t;, + pt^.t;,p + p{p, r^.}r;, + r;,{p, v}p ) 

^[5Pr2] = ( ^-ST^^T^^P + ^PT^,T^,5 - \t^u{S, P}T^, + ^P{S, T^.}T^, 

^[F^T^] = ^( {[iv„ [iv,, T,a]],T;j + ^{[iv^, [iv,, t;j], t,^} 

A[A T ] = i{ —TapAij^A^T*^ + -^AfjToiijT^ijA^^ + -T^Q,[yl^, yljyjT*^, 

4- i 4 ('T T* — T T* ^4 — — T /4 T* 14 

I 1 . 4 . 

4 . 1 . 1 . 

1 2 

+ ^{{-5, t;j, [iP^, Wy]} + t;,], [iP., s]] ) 

A[SA^T] = ~2 ( T*^jA^AyS + (S'A^Ai/T*^^ + A^jS*, T*j^}A^ — 3T*^,{S, Ay}Ai^ 
+ SARIS', Ay\T^j^ — T^y}Ai, + Ay{Ai^, T^^^j^" ) 
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A[PV^T] = i{ -{{T^.,P}, \iV^,W,]} - [[iV^,T^,], [iV,,P]] ) 

A[PA''T] = i{ ^T^^A^A^P + ^PAf^A^Tf,,, - A^{P, T^.}A, + P]A, 

~ A^[A^j P]T^i, — —P[Ai_i,Tfj_p,]A^-\- —A^[A^^T^^]P ) 
A[SVAT] = {-\iV^,S]A,T^, + T^,A,\iV^,S] - \iV^,A,]ST^, + Tf,,S\iV^,A,] 
+ ^[iV^, S]T^,A, - ^A,T^,[iV^, S] - ^[iV^,A,]T^,S 

+ ^ST^,\iV^,A,] + ^[iV^,T^,]SA, - ^A,S[iV^,T^,] 

+ ^\iV^,T^,]A,S -^SA,\iV^,T^,] + '^T^,\iV^,S]A, 

- ^A,[iV^, S]T^, + ^A,[iD^, T^,]S - ^S[iV^, T^,]A, 
+ ^T^.[iP^, A,]S - '^S[iV^, A^]T^, ) 

A[PVAT] = i{ [iv^, p]a,t;, - T;,A,[iv^, p] + [iv^, p]t;,a, 

- A,T;,\iv^, p] - ^[zP^, a,]t;,p + ^pt;,[zp^, a,] 

- \iv^, a,]pt;, + t;,p\iv^, a,] + [iV,,t;^]pa, 

- A,P[iV^,T;^] - ^\iV^,T;j\A,P+^-PA,\iV^,T;j\ 

+ 2t;,\iv^, p]a, - 2A4iv^, p]t;, + '^p\iv^, a,]t;, 

- '^T;,\iV^,A,]P+jp\iV^,T;,]A, - ^A,\iV^,T;,]P ) 



Appendix D. Explicit form of the counterterms in four dimensions 

For all the terms which do not contain iV^ nor A^, one can find a compact expression 
for their counterterms, namely 



W[S, P, T] = -^(M^ + M^^MV + 2MVM7^> (D.l) 



where here M(a;) = S + i^^P + ^iu^vT^v and 5a^{x) = —{ia{x)^5, M(a;)}. The rest of 
the counterterms are given by 
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1 4 2 

+ ...]= 4(1 [zD^, 5][zl?^, 5] + ^[zD^, P][zl?^, P] 
+ ^^A^A^ - PA^PAf, - i\iVf,, S]{P, A^}) 
W[V^T^ + •••]= 4( - ^[A^, T^a] [^., + ^A^A^Taf3Taf3 

- ^[zV„A,]T^pT:p + ^[iV^,T^^]{A^,T:J) 
W[SV^T +•••] = -4i{ - ^{S, T^M^A, + pA^T^^A, + ^i\iV^, W,]{P, T^,} 
+ l[iD^,5][A.,T;j - '^[iV^,A,][S,T;,]-i[iV^,T^,][P,A,]) 
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Table captions 

Table 4.1 Several useful particular values for the coefficients involved in the inverse mass expan- 
sion (4.2) and defined in Eqs. (4.5) and (4.6). The coefficients for the inverse mass 
expanded effective action Eq. (4.14) and the effective current Eq. (4.18) and (4.19) 
are also given. In the odd-dimensional case, the coefficients C,Ki oc^ , Pk^ (^k 
have to be multiplied with an extra ^/tt factor. 

Table 4.2 Explicit formulas for the effective action inverse mass coefficients, and (3^ (see 
formula (4.14)) in terms of (k as given by Eqs. (4.5) and (4.6). 
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Table 4.1 





K<0 i^ = 


K > 




r(-K+l)CKEn=l^ 


{-1)^T{K)Ck 




r(-j^+i)CK Co 






Table 4.2 
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